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New technologies of fiber manufacture and a demand for unusual fiber 
qualities in communication systems have intensified the interest in a 
comprehensive theory of multimode fibers with nonuniform index distribu- 
tions. This paper deals with a general class of circular symmetric profiles 
which comprise the parabolic distribution and the abrupt core-cladding 
index step as special cases. We obtain general results of useful simplicity 
for the impulse response, the mode volume, and the near- and far-field 
power distributions. We suggest a modified parabolic distribution for best 
equalization of mode delay differences. The effective width of the resulting 
impulse is more than four times smaller than that produced by the parabolic 
profile. Of course, practical manufacturing tolerances are likely to in- 
fluence this distribution. A relation is derived between the maximum index 
error and the impulse response. 

I. INTRODUCTION 

Conventional optical fibers consist of a high-index core surrounded 
by a cladding of lower index. The index step between core and cladding 
contains the light inside the core and isolates it from the outer fiber 
surface, whose quality is usually difficult to control. In a more general 
way, inside guidance can be accomplished by any index profile which 
decreases from a maximum inside the fiber to a lower (cladding) value. 
The specific shape of the profile has an effect on the distribution of the 
guided optica] power in the fiber and on the overall loss encountered, 
but, more importantly, the profile profoundly influences the velocities 
of the various propagating modes. A good example is the parabolic 
index distribution which was predicted to nearly equalize the group 
velocities of the propagating modes. 1 ' 2 The Selfoc fiber which closely 
approximates these conditions has indeed since exhibited an extremely 
narrow impulse response. 34 

These effects greatly enhance the chances of multimode fibers to be 
used in optical communication systems. On the other hand, a theory 
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of the interrelations between index profile, impulse response, and 
power distribution is presently only available for the two special cases 
of the uniform and the parabolic core index. This paper provides a 
more general theory and studies a broad class of index profiles po- 
tentially useful in communication applications. The uniform and the 
parabolic profile are special cases within this class. 

Our concern with multimode fibers for communication applications 
allows us to make four simplifying assumptions : 

(i) The index profile is circular symmetric. 
(ii) The core diameter measures hundred wavelengths or more and, 

hence, a great number of modes can propagate. 
(in) The total index change within the guiding core region is only 

a few hundredths, so the propagating modes can be considered 

essentially as transverse electromagnetic. 6 
(iv) Index variations within the distance of a wavelength are 

negligible, and the conditions of geometrical optics (or the 

zeroth order of the WKB method) apply. 

Except for these four restrictions and the requirement of guidance, 
the index profile can be of the most general form. It can, for example, 
have an index depression in the center and one or several ring-shaped 
index maxima. 6 

For the sake of clarity, this paper is restricted to the simpler type 
of profile illustrated in Fig. 1. We assume the index profile will decrease 
monotonically from the center and converge into a flat cladding region 
which guarantees isolation from the outside surface. 
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Fig. 1— Cross-sectional sketch of circular symmetric index profile in multimode 
fiber. 
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Fig- 2 — Wave vector diagram in the propagating region of a multimode fiber. 

Apart from the index profile, there are, of course, other influences 
which affect the impulse response and the optical power distribution 
inside the fiber. Mode excitation, loss differences in the process of 
propagation, and coupling among the modes play a part. To isolate 
the effect of the index profile, we assume here the ideal case of uniform 
loss, absence of coupling, and equal and simultaneous excitation of all 
propagating modes at the input. For the computation of the impulse 
response, the input is assumed to be an infinitely narrow pulse of unit 
energy. 



II. MODE DESIGNATION AND MODE COUNT 

All guided modes are essentially transverse electromagnetic and, 
with some proviso, can be decomposed into linearly polarized pairs. 5 - 7 
Because of the circular symmetry of the index n, the modes have a 
circular periodicity and can be identified in the conventional way by 
an azimuthal order number v. To characterize the radial field distribu- 
tion, we need an additional mode number m- The propagation constant 
of a particular mode (/*, v) can then be approximately determined 
by the WKB method. 68 Figures 2 and 3 give a physical description 
of these relationships. In Fig. 2, the local wave number 



k(r) = 2irn(r)/X 



(1) 
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Fig. 3 — Sketch defining regions of periodic and aperiodic field characteristics of a 
mode of azimuthal order v. 

is decomposed into its components in a cylindrical coordinate system 
(r, <f>, z). The unknown radial component becomes 

u(r) = [fc 2 (r) - 2 - Wf. (2) 

Given j8 and v, we can find two radii Ri and R 2 , at which u(r) vanishes 
(see Fig. 3). These radii define a ring-shaped region within which u is 
real, causing a radial periodicity of the mode field. Outside this region, 
the field is aperiodic. 

Radially decreasing (or evanescent) field conditions obtain outside, 
when the phase inside (approximately) adds up to an integer number 
of half periods between Ri and Ri. Consequently, if y. designates this 
number of half periods, 



M7r = (*' U (r)dr - P [/c 2 (r) - (P - v 2 /r^dr. 

J R l J R x 



(3) 



We would have obtained the same result by way of the WKB method, 
with the only difference that ft and v 2 would be replaced by n + | and 
i> 2 + I. These corrections are important in the case of small n or v, 
and particularly for the fundamental mode which has n = v = 0. On 
the other hand, to obtain a general view of the mode structure, we 
can ignore the £-terms as long as we refrain from discussing individual 
low-order modes. 

For the purpose of a total mode count, let us consider the limits of 
H t v, and 0. The requirement of evanescent field conditions in the 
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cladding (index n c in Fig. 1) limit /3 to a minimum value 

C = 2*n e /\. (4) 

Modes with smaller /3 find propagating conditions in the cladding and 
are no longer bounded by the core profile. Condition (4) defines mode 
cutoff. The largest value for v results for = /3 C and n = 0, and 
alternatively n is largest for = ft e and v = 0. We obtain the total 
number of modes M from a summation of (3) over all v from to 
"max- If v mKX is a large number, we may consider v a continuous variable 
and replace the sum by an integral. In this case, 

M = - / [fc 2 (r) - /3 2 - vyr^drdv. (5) 

7T JO J Ri(v) 

The factor 4 in front of the expression allows for the fact that each 
combination /x, v designates a (degenerate) group of four modes of 
different polarization or orientation. 5 Figure 4 illustrates the area of 
the double integration indicated in (5). A change of order in the 
integration leads to 

= - / / (*"-#- */t*)*dvdr, (6) 

IT JO JO 



M 



where a is the radius at which the index n(r) reaches the cladding 
value n e . Integrating (6) with respect to v yields 



M 



= f* LHr) - 0Sfrdr = ( y)' f* 2 (r) - n\^rdr . (7) 

For small index differences, the integral represents the volume under 
the (circular symmetric) profile plot. It may be worth noting, though, 
that the substance of this relation is not limited to circular symmetry. 




Ri(i') R 2 („] 

RADIUS »- 

Fig. 4 — Region of double integration in eq. 5. 
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RADIUS r *- 

Fig, 5_a few of the index profiles defined by n = n [l - 2A(r/a) a ]* for small A. 



For later use we write (7) in the somewhat different form 

•fl2(0) 



/•fl2(0) 

m(B) = / [fc 2 (r) - pyrdr, 
Jo 



(8) 



where m(j3) denotes the number of modes having a propagation con- 
stant larger than 0. The upper limit R 2 (0) of the integration is the 
radius at which k(r) = /?. 

Let us now consider a particular class of profiles defined by 



n(r) = 



n [\ - 2A(r/a)°y for r<a 



(9) 



n [\ - 2A]* 



for r > a , 



where a is a parameter between 1 and « . Figure 5 illustrates the cases 
a = 1, 2, 4, 10, and °o. All profiles reach a constant cladding value at 
r = a. The core profile has a cone shape for a = 1, becomes nearly 
parabolic for a = 2, and converges to the case of the step profile for 
a = oo. Using (1) we introduce (9) into (8) and obtain 



/// 



tf ).^-A_ 2 («^) 



|2 \ (2/a)+l 



(10) 
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where 

k = 2irn /\. (11) 

For fi = p e from (4), the total mode number becomes 

M = -^aW A. (12) 

It is proportional to the index difference and the core cross section. 
The uniform profile accepts twice as many modes as the parabolic one 
and three times more than the cone-shaped one. 

III. IMPULSE RESPONSE 

Consider all modes to be excited by the same narrow pulse at the 
input. Each mode transports an equal amount of energy to the fiber 
end. The individual pulses are expected to suffer a certain distortion, 
depending on the 0-co characteristic of each mode and dispersion in the 
dielectric. We assume, however, that the resultant broadening is small, 
or at least not much larger than the group delay differences between 
adjacent modes. Because of this effect and other limitations in the 
system response, the pulses from individual modes are likely to fuse 
into one continuous output pulse called the impulse response. Since 
all modes carry the same energy, the power profile of the impulse 
response is equal to the mode density per unit time interval. In the 
following theory, the continuity of the impulse response results not 
from the broadening of the individual mode responses, but from the 
assumption that n and v are continuous functions. 

The straightforward method of computing the impulse response 
starts from (3) to find the propagation constant /3 for each pair, n, v. 
The group delay in a fiber of length L is then 

«,')-^, (13) 

where c is the vacuum velocity of light. A simplification of this ap- 
proach for the purpose of numerical computations is indicated in the 
appendix. Once t(h, v) is known, the impulse response results from 
a count of the combinations n, v which arrive between t and t + dr. 
This number plotted versus t then constitutes the impulse response. 

For the special class of profiles defined by (9), group delay and 
impulse response can be computed in a much simpler way. First we 
postulate that, in this case, the relation between r and /3 according to 
(13) is independent of n and v. If this holds — and we shall prove it 
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later with the help of eq. (16) — we can replace by r in (3) and still 
perform the same integration over v which led to (8) and, more specifi- 
cally, to (10). Solving the result of this integration for t yields 

r - £2! < \ K - ( 2 4^)"" + "(2 4 ff«.«l'. (14) 
c dk l \ a? a / J 

This result can easily be verified by solving (10) for and introducing 
it into (13). With the help of (10) and the abbreviation 

8 - 1(1 - 0Vfc» o ) , (15) 
eq. (14) takes tho form 

Ln l -4V(« + 2) , lfi v 

T = ^ (1 - 25)* (lb) 

This expression proves indeed to depend on alone (and not explicitly 
on w), thus justifying the approach chosen. 

To obtain the impulse response, we can now introduce (16) into 
(10) and differentiate with respect to r. Although this is not difficult 
to do, it leads to rather unwieldy expressions. We shall therefore merely 
consider some special cases of interest. To normalize the impulse 
response for total unit energy, we divide (10) by (12) and obtain 

* /iy»'-*\ (n) 

M W 

Furthermore, since 5 can at most assume the value A (for /3 = /3 C ) 
and is therefore small compared to unity within the scope of our theory, 
we develop (16) into a power series in terms of 5 and obtain 



Ln / a - 2 3a - 2 

' = — ( 1 + ^r+2 5 + 7r+2 



|). (18) 



We relate t to the total propagation time Ln /c and introduce a new 
time reference, which ignores the delay common to all modes. Hence, 



re a - 2 . . 3a - 2 5 2 



,. " l.«JZ-« + ~-i. (19) 

Ln a. + 2 a + 2 Z 

In this time frame, the fundamental mode arrives at t = 0. 

As long as a is not too close to 2, the linear term in (19) dominates. 
Therefore, 



5 = 



" 2 t except for a«2 



a - 2 

V2J for a = 2. 



(20) 
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Insert this into (17) and differentiate with respect to t to obtain the 
impulse response 



1 dm 
M~dl 



^ a + 2 \ a + 2 1 

a | a — 2 A 

2^ 

A 2 



(2/a)+l 



\t\ 2 "* except for a«2 



(21) 



for 



= 2. 



As 6 varies from to A, the time t changes from to 

fa - 2 

except for a tt 2 



T = 



a: + 2 

A* 

2 



(22) 



for a = 2. 



Outside of this time interval, the impulse response is zero. Figure 6 
shows plots of (21) for the profiles sketched in Fig. 5. A change from 
a = * to a = 10, which implies a relatively small change in the 
profile, narrows the impulse response by £. The response becomes 
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Fig. 6 — Impulse response of multimode fibers having the profiles of Fig. 5. 
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extremely narrow for a « 2, then broadens again, as a decreases 
further. For a < 2, the high-order modes overtake the fundamental 
and arrive earlier. 

In the vicinity of a = 2, where both terms of (19) contribute, the 
impulse response is a rather complicated function. The most interest- 
ing of these cases is the one for which the impulse response has the 
narrowest possible width. This optimum condition arises for 

a opt = 2 - 2A , (23) 

which yields 

t = J(* - AS). (24) 

In this case, the modes of highest and lowest order both arrive at the 
same time t = 0; all other modes are faster, the fastest one being 
determined by 5 = A/2. It arrives at 

.--£ (25) 

Equation (24) has two solutions for 8. Hence, (17) yields two values 
for the same t, indicating that two mode groups, a high and a lower 
order, contribute to the impulse response at every particular instant 
in time. By introducing 8 into (17), differentiating with respect to t, 
and then adding the two contributions, we find the impulse response 



A 2 V ^ AV 



"' (26) 



This function is plotted in Fig. 7. It peaks at t = - A 2 /8 and decreases 
towards t = 0. Because of the normalization introduced in (19), the 
absolute temporal width is 

Ln A\ (27) 

c 8 

The time slot in which a pulse of this kind can be transmitted is 
narrower than that, because 70 percent of the power is concentrated 
in the first half of the interval (27). 

A practical implementation must, of course, allow for a certain 
tolerance or error in the profile, as a result of which the total width 
of the impulse response is likely to exceed (27). To obtain some in- 
dication of the pulse broadening as a result of this index deviation, we 
assume that the erroneous profile is still of the type (9), but has 

a = a pt + da. (28) 
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Fig. 7 — Impulse response in the case of optimal profile shape. 



The maximum index deviation from the optimum profile then appears 
approximately at 

r = ae-l (29) 

and has the value 

n A 



dn m&lc = di 



2e 



(30) 



where e is the base of the natural logarithm. As a result of this profile 
error, the normalized width of the impulse response becomes 



or, in absolute terms, 



|(A + i[da|) 2 



Ln, 
~8c 



A + 



n„A 



[dn, 



(31) 



(32) 



Consider a guide with a maximum index n = 1.5 and an index 
variation A = 2 percent. If the profile is optimal, mode delay should 
produce an effective broadening of only 0.25 ns/km. An index devia- 
tion of 10 -4 from the optimal profile increases the broadening to 
0.53 ns/km. 

IV. NEAR- AND FAR-FIELD POWER DISTRIBUTION 

We take again into account the fact that the core cross section mea- 
sures many wavelengths in diameter. If this cross section is illuminated 
by an incoherent source (exciting all modes uniformly), the power 
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incident per unit solid angle at any point in the cross section is con- 
stant. To compute the power accepted by the fiber, we merely have to 
know the solid angle of acceptance at any point. We find this angle 
from the wave vector diagram of Fig. 2, which yields 



cos 6(r) = 0/k(r). 
The maximum angle d c results for = fi e ; hence, 



cos c (r) = 



k(r) n(r) 



(33) 



(34) 



Using this relation, we can define a local numerical aperture at the fiber 
front face 

A(r) = n(r) sin 9 e (r) = [n^r) - n 2 c ]*. (35) 



The power accepted at r is then 



. ... A*(r) . n*(r) - n\ 

P(r) = p(0) 1 ™=p(0) n2(0) _ n2; 



(36) 



If all modes propagate equally attenuated and without coupling, the 




RADIUS r +- 

Fig. 8 — Power distribution in the core of multimode fibers having the profiles of 
Fig. 5. 
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Fig. 9 — Power distribution in the far field of multimode fibers having the profiles 
of Fig. 5. 

same power distribution should hold for the fiber end face. The class 
of profiles described by (9) has 

A(r) = n (2A)*[l - (r/a)«y (37) 

and 

p(r) = p(0)[l - (r/a)°J (38) 

The agreement between the profile plots (Fig. 5) and the near-field 
power plots (Fig. 8) is not a coincidence, but holds in general as long 
as the total index variation is small. 

Under the conditions assumed here, every incremental area of the 
core cross section at the fiber end uniformly illuminates its cone of 
acceptance. For this reason, all those areas that have a numerical 
aperture 

A (r) ^ sin 6 (39) 

contribute equally to the far-field power at 6. For the class of profiles 
described by (9), the areas contributing to 6 are within a circle whose 
radius is obtained by solving (37) for r. Consequently, 

ro - p«> (i - gy (40) 
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is the far-field power distribution. Figure 9 shows a plot of (40) for 
the profiles of Fig. 5. The uniform illumination for a = <» changes 
to a parabolic distribution for a = 2. All plots must be understood as 
the average power expected under the idealized conditions mentioned 
earlier. Monochromatic mode excitation results in mode interference 
phenomena and a local fine structure, which can greatly modify the 
average distribution considered here. 

V. CONCLUSIONS 

By assuming somewhat idealized conditions for mode excitation, 
coupling, and loss in a multimode fiber, we can isolate the influence 
of the index profile upon mode volume, near- and far-field power dis- 
tribution, group delay, and impulse response. Surprisingly simple 
relations exist for a special class of profiles which comprises most 
multimode fibers of interest. Particular attention is given to a near- 
parabolic profile which accomplishes optimal delay equalization of all 
modes. If the (relative) index difference between center and periphery 
of this profile is A, mode delay broadens the impulse response by a 
fraction A 2 /8 of the total propagation time. This amounts to about 
0.25 ns/km for A = 2 percent. On the other hand, an index deviation of 
10 -4 from the optimal profile increases the broadening to 0.53 ns/km. 

APPENDIX 

Some Further Relations for the Group Velocity 

The numerical evaluation of /? as a function of n, v, and k and its 
subsequent differentiation to obtain r are usually tedious and time- 
consuming. A substantial simplification results from a direct computa- 
tion of r by applying the operation 

Ln dfi/dko 



to (3). The result is 



c dp/ dp 

I k(r)n(r)dr/u(r) 

jRi 



(41) 



'-v^ (i2) 



fiil 

/ (3dr/u(r) 

J Ri 



To understand the physical significance of this relation, consider a 
ray propagating along the fiber core of Fig. 2 in such a way that it has 
the direction of k (r) at r. A line element along this ray is 

ds = (dr* + rHtf 4- dz 2 )* (43) 
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and therefore 

ds k(r) dz B 

dr = W) ^ dr = WY (44) 

The condition u(r) = dr = at R x and R 2 indicates a reflection (turn- 
around) of the ray. The ray performs periodic undulations between 
Ri and R 2 , simultaneously moving sideways in a helical fashion. By 
introducing (44) into (42), we obtain 

<f> n(r)ds/c 
-L^- . t (45) 



/ 



dz 



where f denotes integration over a full period of the ray. The de- 
nominator describes the axial length of one ray period, and the nu- 
merator the propagation time along the ray within this length. 
Multiplied by the fiber length, this ratio yields the total group delay. 
This result emphasizes the equivalence between ray theory and the 
zeroth-order WKB approach followed in this paper. 

Within this order of approximation, the only quantities that depend 
on the wavelength are the mode numbers. Normalization of these 
numbers and subsequent transition to continuous variables eliminates 
the wavelength entirely ; group velocity and impulse response are then 
independent of wavelength. More specifically, if we write 

p - v./ak and a = v/dk (46) 

and 

n = n„{\ - 2d(r)y , (47) 

eq. (3) assumes the form 

and (42) becomes 

t„ f (1 - 2d)dr/[28 - 2d - («ro/r) 8 ]* 

f ' dr/[25 - 2d - (<m/r) 2 ]* 



-' (1 ~ 2«)"» :LJ ^— F W (49) 



U 



These two equations are sufficient to calculate group velocity and 
impulse response in the case of large mode numbers. 
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